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ABSTRACT 

An adaptive oontrol system operating on submarine depth rate 
with constraints on the maximum values of stern plane deflection 
angle and hull pitch angle is formulated using a linear compensation 
technique. The system is adaptive in the sense that it is designed 
to give near uniform dynamic response over ranges of values of various 
hull geometry and speed parameters. The formulated system is simulated 
on a Reeves Electronic Analog Computer, various adjustments made and 
justified, and results shown in plots of stern plane angle, pitch angle 
and depth rate error. A simple system for predicting pull-out time is 
suggested which, when associated with the proposed system, would 
facilitate rapid, well-damped depth changes. 
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Symbols and Nomenclature 



The nomenclature defined in T & R Bulletin No. 1-5 of the Society 
of Naval Architects and Marine Engineers is used herein where applic- 
able. A right-handed cartesian coordinate system is employed with the 
origin at the center of gravity of the hull, the positive x axis dir- 
ected forward, the positive Z axis directed down normal to the base 
line, and angles measured positive in the counter-clockwise direction. 
Pitch angle and depth are measured relative to a similar system fixed 
in space and oriented with the positive Z axis directed vertically 
down. Quantities measured relative to this fixed system are designated 
by the subscript g. All quantities are non-dimensionalized, as indi- 
cated by the prime, by dividing by appropriate powers of the length of 
the hull and linear velocity, the mass and mass density. 

non-dimensionalizing 



symbol 

I 

y 


I* 

y 


formula 

I 

y 


definition 

Moment of inertia of 
hull about Y axis 


k 

y 


k' 

y 


rs 


k 

y 

L 


Radius of gyration of 
hull and added mass of 
hull about T axis 




k ' 
*"2 


SB 


^2 


Lateral added mass 
coefficient 








L 


k' 








Added moment of inertia 
coefficient 


L 


L' 


=: 


1 


Characteristic length 
of hull 


m 


m' 


rs 


m 

•> 


Mass of hull 













M 

q 



M. 

q 



M 

w 



M. 

w 




q 

4 



u 



w 






z. 

q 



z 



w 



M' 

q 



M 



i/oL^U 



M, ' 

4 



M. 



hpl' 



M' 

w 



M 



w 



^-/oL^ 



Derivative of moment 
component with respect 
to angular velocity 
component q 

Derivative of moment 
component with respect 
to angular acceleration 
component 4 

Derivative of moment 
component relative to 
normal velocity component 
w 



M. • 

w 




w 

%s 



u^ 



Derivative of moment 
component relative to 
normal acceleration 

Derivative of moment 
component relative to 
stern plane angle 




Angular velocity com- 
ponent relative to y axis 

Angular acceleration 
component relative to 
y axis 



u’ 



w’ 



w* 



"4' 



w 

IT 



i/OL^ 



Z’ 

W 



w 



i/OL u 



Velocity of origin of 
body axes relative to 
fixed axes in feet per 
second 

Component of u along Z 
axis of body 

Component of acceleration 
of origin of body axes 
along Z axis 

Derivative of normal force 
component with respect to 
angular acceleration 
component 

Derivative of normal force 
component with respect to 
normal velocity component 
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%s 


Derivative of normal force 


i/oLV 


component with respect to 
stern plane angle 




metacentric height 


L 






stern plane angle 


e rad 


Pitch angle as measured 
between X axes of body 
and fixed coordinate 
systems 


1 


mass density of water 




time in seconds 


subscripts 





desired or commanded quantity 



actual or output quantity 
relative to fixed coordinate system 



I 



INTRODUCTION 

A. General Aspects of the Depth Control Problem 

Prior to the closing years of World War II, submarine design features 
were based primarily on surface performance. As a direct result of this 
philosophy, submerged speeds, and hence maneuverability, were slow and 
sluggish. The top speed of a typical boat of this period was only nine 
knots at the half-hour battery rate. For this reason the crude two- 
plane sman depth control system generally employed since the very earliest 
days of submarine construction was entirely adequate. 

The improved methods of anti-submarine warfare brought about by 
developments in sonar and radar, the use of carrier-bourne aircraft, 
and the construction of fast anti-submarine escort ships made higher 
submerged speeds imperative. Belatedly, submerged performance became 
the primary consideration. By greatly reducing the amount of super- 
structure, and fairing the bridge and periscope sheers, hull resistance 
submerged was reduced by almost a factor of two. This, coupled with a 
doubling of battery capacity, increased submerged speed by a factor of 
two. These higher speeds indicated clearly the need for faster, more 
accurate control of depth. This need is easily recognized when one 
considers that submarines of this period would collapse at depths little 
in excess of one ship length. Even present-day submarines with their 
much higher speeds must confine their motions to a layer of water only 
a few ship-lengths deep. 

The advent of nuclear propulsion plants and the low resistance. 
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“body of revolution hull form developed in the last decade have greatly 
magnified this depth control problem. Let us consider the situation con- 
fronting the submarine pilot controlling the stern planes when only stern 
planes are used to control pitch angle and depth. This situation is 
‘adequately represented by the following block diagram. 




Zgo actual depth 0 pitch angle 

Pig. I 

Prom this diagram it is evident that there are three integrations be- 
tween the immediate e^^ect of the operator's action and the desired 
result. If the operator is to close the control loop in an optimum 
manner, he must anticipate the results of these integrations - a matter 
difficult in itself, but greatly complicated by the long lag between 0 

and 5 introduced by the la*ge hydrodynamic mass and small hydrodynamic 
s 

stiffness of the hull. The operator attempts to overcome this diffi- 
culty by introducing a large initial § . This gives him a more rapid 

s 

indication of the result of his control action, but due to the lag 
mentioned above, he has no timely indication of when to remove or re- 
verse plane angle. This results in overshooting the desired objective 
- the behavior of a poorly damped system. Thus it is evident that the 
operator is required to perform two functions simultaneously in closing 
the control loop: he must speed up the system response, and he must 
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provide adequate damping characteristics. This latter skill requires 
a certain amount of prescience highly dependent not only on the skill 
of the operator, but also on his physical condition. In this sense, 
operator reaction is a random variable introduced into the control i>ro- 
blem which precludes absolutely consistent performance. 

Two basic approaches have been followed to alleviate or avoid these 
difficulties in depth control. The first, and most apparent, is replace- 
ment of the human operator in the control loop by an automatic plane 
angle computer. The second approach, and the one followed in this in- 
vestigation, involves relocation of the human operator outside of the 
closed-loop portion of the system, 

B, Development of Automatic Depth Control Systems 

The German Navy was the first to practice the new philosophy of 
submarine design. Their type XXI boats constructed just before the end 
of World War II in Europe were capable of about twenty knots submerged. 
These ships clearly indicated the desirability of an automatic depth 
control system. An attempt was made to meet this need. The system 
which evolved was really two very similar sub-systems - one for each 

5 

set of planes. Plane angles were computed as functions of speed, 

depth, and pitch angle. Depth and pitch angle signals were derived 
from pressure gauges and pendulums respectively. Variable gain settings 
in each of the sub-systems permitted ratioing bow and stern plane angles 
as a function of speed and depth-keeping accuracy requirements. The 
controlled quantity was depth, so that the overall system was basically 
a depth regulator rather than an autopilot capable of radical maneuver- 
ing. 
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The United States, profiting from German experience, began modifying 

the top-hamper and battery capacity of their SS-475 class fleet boats in 

* 

the late 1940' s. Coincident with this Guppy conversion, development of 
a semi-automatic depth-keeping system was initiated. The first such 
system was built by Askania Regulator Company and installed on the 
U.S.S. Irex (SS-482) in 1948. Although very similar to the German 
system in most respects, it controlled not only depth, but also boat 
angle or pitch angle. With this added feature, it provided more man- 

5 

euvering capability than the German system. The first of these systems 
were semi-automatic in that the output was an indicated plane angle 
which the planesman was required to match. Subsequent systems were pro- 
vided with an automatic follow-up device making them fully automatic. 

With the advent of Albacore and submerged speeds in excess of 
twenty knots, these early systems became inadequate. Present systems 
are fully automatic, and similar to aircraft autopilots. The Sperry 
System, installed on most new construction, determines stern and sail 
plane angles from signals in depth, depth rate, pitch angle, pitch rate, 
and pitch acceleration to keep the boat on a digitally computed optimum 
trajectory. 

C. Developments in Manual Control Systems 

The high submerged speeds and limited depth range of modern sub- 
marines make the poor coordination and large time lags inherent in the 

*The Guppy (Greater Underwater Propulsive Power) conversion involved 
streamlining the superstructure and doubling the battery capacity 
of the SS-475 class boats. These two modifications resulted in a 
200^ increase in submerged speed. The resulting ships bear an un- 
canny resemblance to the type XXI boats. 
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old two-operator, manual control system intolerable. The first step 
away from the old system was incorporated in Albacore with the substi- 
tution of an aircraft-type control column for the old wheel type de- 
vices. Two sticks were provided. The main unit is a conventional 
aircraft steering column in which wheel rotation controls rudder posi- 
tion, and fore and aft motion of the entire column controls bow plane 
angle. The second stick controls stern plane angle. In the principal 
mode of operation, stern plane angle command is derived from the bow 
plane signal via a variable gain amplifier, thus providing single- 
operator ratioed control. This is the mode of operation greatly pre- 
ferred by the operators. 

This basic system has proven satisfactory in most respects. Sub- 
sequent developments have been primarily concerned with providing the 
operator with the required feed-back information - depth error and boat 
angle - in the most efficient manner. The first major improvement in 
this line was the Combined Instrument Panel which displays all the 

information required by the operator in a cluster of dials arranged 

8 11 

in an optimum manner from the human engineering point of view. ^ A 
second, and rather novel, development is the Contact Analog Display 

g 

developed by the Electric Boat Division of General Dynamics. In es- 
sence, this is a television screen picture of the path the operator is 
required to follow. A three dimensional effect is provided by a fic- 
ticious horizon and narrowing-with-distance perspective of the pic- 
tured road. This is the only indicator the pilot must concentrate his 
attention upon. His task is simply to keep the white line of the road 
properly positioned on the screen. 

In anticipation of the performance capabilities to be required of 
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future high speed submarines, the Bureau of Ships has initiated an ex- 
tensive study of the over-all submarine control problem. The title of 
this study is the Submarine Integrated Control Program (SUBIC). The 
Contact Analog mentioned above originated in this program. As a part 
of this investigation, Cornell Aeronautical Laboratories conducted a 
study of the data flow and data processing requirements in effective 

7 

submarine control. One of the major results of this study was an 
indication of the desirability of locating the human operator external 
to a multiple closed-loop system in which the outer loop closure is on 
depth rate. The advantages of this scheme are readily apparent. In 
addition to minimizing the effects of erratic human behavior and sim- 
plifying the operators' problem by avoiding two of the integrations 
mentioned in (a) above, the system is easily made fully automatic by 
substituting a digital computer or programmer in place of the human 
operator. This system is represented by the following block diagram. 
The K's indicated are in general non-linear functions of ship speed and 




Pig. II 

the geometrical parameters of the hull. 
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D. The System to be Investigated 

The major difficulty in submarine control system design is the lack 
of precise information concerning the values of the various hydroframe 
parameters. At the present time, there is no satisfactory method of es- 
timating these parameter values within acceptable accuracy from the 

geometry of the hull as delineated on the lines drawing. Recourse to 

12 

model and full-scale tests is unavoidable. Further, in order to 
obtain acceptable performance over the speed range of a given hull, the 
characteristics of a control system such as that shown in Fig. II must 
be varied in some pre-determined, generally non-linear manner with 
speed. It is with this problem that this thesis is concerned. 

The basic system under consideration is shown in Fig. III. It is 
evident from this diagram that no provision is made for a steady state 




Fig. Ill 

plane angle in a steady dive. This feature would require an auxiliary 
forward loop from the input to the plane angle summer. This matter is 
ignored in the following investigation because the effect of such a for- 
ward loop is small and decreases with speed roughly as the reciprocal of 
the speed cubed. The effects of this simplification will be a slight 
increase in steady state depth rate error at slow speeds. Also, only 
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stern plane control is considered. 

The objective of this investigation is to formulate suitable com- 
pensation functions K ^ and K G ^ that will yield a specified per- 

CX CX Gc. G d. 

formance of the closed loop system over a specified speed range of a 
given hull configuration, and at the same time minimize the effects of 
deviations of certain important hull parameters from their specified 
values. In this sense, the proposed system is an adaptive system. 

The problem will be treated on a non-dimensional basis using the 
linearized force and moment equations of a submerged body in the ver- 
tical plane, and standard nomenclature as specified in reference 
The system will be simulated on a REAC analog computer to determine its 
transient and steady state response to a step input of depth rate under 
various conditions. 
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PROCEDURE 

A. Scope 

This investigation is concerned with the formulation of a suitable 
closed-loop system for controlling the depth rate of a deeply submerged 
submarine subject to constraints on the maximum values of hull pitch 
angle and plane deflection angle. The principal objective is a system 
whose dynamic response to an ordered change of depth rate is relatively 
independent of variations in ship speed and hull geometry. A linear 
design technique will be employed, and the performance of the resulting 
system analyzed on an electronic analog computer. 

The basic idea underlying the approach to be followed is that 
hydroframe parameters may be considered to fluctuate from hull to hull 
as well as for a given hull about certain normal values. Similarly, 
ship speed may be considered as a deviation from some normal speed. 

The effects of these deviations is a change in the resonant frequency 
of the hull and hence of the entire control system. If when its para- 
meters have their so-called normal values, the overall system is so 
designed that its resonant frequency is much higher than any frequency 
to which it is normally required to respond, then the change in natural 
frequency brought about by parameter variations will have a small effect 
on system performance. 

The design specifications for the proposed system are as follows: 

1. a non-dimensional natural frequency of the **normal** system 
in the vicinity of 5 with a damping ratio of about .75* 
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2. no more than 15 ^ overshoot in pitch angle.* 

3. no more than 10 "^ overshoot in depth rate. 

4. a rise time to 10° pitch angle for large changes of depth 
rate close to 1 second in non-dimensional time. 

The non-dimensional form of the equations of motion of the hull 
are employed since this facilitates separation of the dominant para- 
meters from those of lesser importance o 
B. Equations of Motion of the Hull^ ^ 

Omitting the cross-acceleration terms M. w' and Z.'., which are 
° w q q' 

2 

usually quite small, the linearized equations of motion of a submerged 
body in the vertical plane may be written:^ 

(1) (m’-Z.') Z’ w' - (m’ + Z’) q' - Z1 § 



w 



w 



5s ° 



(2) -M' w' + (!' - M' ) 4' - M' q' - Z1 e = M« g 
'■'^w 'y q q 6 s 

s 



The coefficients of the acceleration terms w' and q' are inter- 
related through the following equations: 

I* » (k* m' 

y y 

Z. ' = -k„m' 
w 2 

M.' c -k’ I’ 

4 y 

k’ is simply the non-dimensional radius of gyration of the hull about 
the y axis, represent the added mass of the hull. They are 

usually taken as the theoretical values applying to the ellipsoid of 
revolution approximated by the hull form.^^ 

*It should be noted that the desired glide angle is implicitly included 
in the commanded depth rate since Z = w - u6. Assuming a steady glide 
with zero normal force, w » o and Z® = - U 0. In non-dimensional nota- 



tion, Z' = - e 

* g 



g 
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Landweber and Johnson of the David Taylor Model Basin have pro- 
posed the following relations among the coefficients of the velocity 

terms for the case under consideration where the only appendage is a 
12 

tail surface. 

M* o m' - X' Z' 

w 1 w 



Z' 

<1 



XI Z' 
2 w 



M' = (X!)^ Z’ 

q ^ 3' w 

These relations, although only approximate, give predictions of the 
indicated parajneters which are reasonably close to measured values. 
xj may be interpreted as the non-dimensional effective moment arm at 
which the lift due to angle of attack acts. The m' in this relation 
is taken as an approximation to the Munk moment which appears in the 
potential flow about the body. X^ may be regarded as the non-dimen- 
sional effective radius through which the angular velocity q is related 
to an equivalent normal velocity, i.e. w' = X^ q' . (X^) may be inter- 

preted as the effective product of the radius relating angular and 
linear velocities and the moment arm relating moment and force due to 
angle of attack. Although , such a product is not necessarily a square, 
it seems nearly so and is so defined for convenience.^*^ 

The coefficient of 6 in equation (2) may be reduced to a speed 
parameter as follows: m| Zg 

Z • 

- ^ 



i/oL"U" L 



8 Z 



B 



ra’ 



rrr) ^ 



u" 



u 



This speed parameter U is then dependent not only on the ship speed U, 
but also on the metacentric height which is a measure of the ship's 
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longitudinal static stability. Incorporation of the minus sign in the 
definition of l/ff is done since is always a negative quantity for 
ships with positive static stability. 

The formulation is further simplified with the introduction of the 
two quantities 



K' 

s 



Z’ 

<^s 



X' 

S 









K’ may be regarded as a measure of the control effectiveness or as the 
s 



ratio of the lift force generated by the stern planes to the lift force 



generated by the fixed surfaces of the hull due to angle of attack. X' 



is the non-dimensional effective radius at which the normal force gen- 
erated by the stern plane acts. The negative sign arises from the 
coordinate system employed. 

Making all of the above substitutions in equations (l) and (2) 
after dividing through by ra' yields the equations of motion of the 
hull with which we shall be concerned. 



(3) (1 + k2) w' - - ( 1 + X* q« = ^ ^ 

m* m’ s 



m’ 



(4) -(1-X«z; )w- + (kp'" (1 + k«) 4’ - (X 3 ‘) q'+9_= -X^i^'z; ^ 






m* 






The utility of this formulation lies in the fact that U is the 
only speed-dependent parameter of importance, and the number of geome- 
try-dependent parameters has been greatly reduced. If it is assumed 

that Z varies as the first power of the speed, which is not unreason- 
w 
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able, then it can be seen that Z^’/m* is independent of speed from the 



following: 



2w- 



w 



u 



Z L 
w 



m 



i/OL- 



Um 



For hulls of normal design, the ranges of variation of k* , and 
k’ -are very limited. Similarly, the effective moment arms and radii 

y 

2 

, ^2* ’ confined to sufficiently limited ranges of 

values that they may be considered to he of secondary importance in 
this investigation. Accordingly, these parameters are assigned the 
following normal values and henceforth treated as constants. 



x^’ = -0.5 


(ky’)^ = 0.045 


X2’ = 0.5 


k^ = 0.933 


(x^’)^ = 0.25 


k’ = 0.805 



X ’ = -0.5 

s 

Then the dynamic characteristics of the hulls to he considered 

, 

depend upon the three primary parameters U, Zw*/®’ » and K *. As 

s 

noted above, the metacentric height Z^ also enters implicitly as a 
forth major parameter. These primary parameters will he considered to 
fluctuate about assigned normal values as follows: 



U 


0 

• 

D 


+ 


4.0 


Zw' 


= .-1.5 


+ 


0.5 


m’ 








K’ 

s 


» 0.3 


-f 


0.1 



It is considered that the indicated ranges of Zw'/m* and K ’ 

s 

include most body of revolution hulls existing c^r to he constructed. 
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The actual speed range indicated by the range of U depends, of course, 
on the value of metacentric height selected. This figure normally lies 
in the range -O .5 to -2.0 feet. For = -0.5 feet, the indicated range 

/sJ 

of U corresponds to a speed range of 2.3 to 21 knots. For = - 2.0 
feet, this range corresponds to speeds from 4.7 to 41 knots. Again, it 
is considered that the ra^ge of U considered includes all combina- 
tions of Z and U of practical significance. 

.D 

C. Formulation of the Proposed System 

By applying Laplace transform operations, in which the operator "s" 
implies differentiation with respect to non-dimensional time , to 
equations ( 3 ) and (4)> one obtains the following transfer functions for 



the hull . 

(5) e( s ) 



(Ej^s + F^) s 
_ _ 

As +Bs + Cs + D 



( 6 ) 






(s) 






Ni s^ - s - N3 
As^ + Bs + Cs + B 



where 



A= (l+k^Kky*)^ d+k’) 



B = -Z 



^ [(k'y)^ U+k') + (l+k^) 

/Z ' \^r . n z . •) + 1 + k 

{^) Vj - ^ — 



- 1 



D - - Z ' 1 

w 

m' -2 



-d+kg) Z 



w 



m'' 



Z ’ K ’ 
1 w s 

m* 



Z ' 
w 



(x^' - x^') + 1 
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Equations (5) and (6) are rewritten in standard from as follows 



m 



K^( s+N) s 



3 2 

s + J^S + J- s + J 
d 1 o 



(8) Z • (s) 

~§TiT 



(s% £^3+ P^) 

3 r 2 ' 7" 

S + J-S + J, s+ J 
2 1 o 



where 



IB E^ 

a“ 

N = P^ 

"e7 



K, 



E, 



B 

T 



P, 




Substituting the assigned values of the secondary parameters 
given above, these coefficients are expressed as functions of the 
primary parameters as follows: 



K, - 6.16 K ' Z ’ 

1 s w 

m' 



= 0.517 K • Z ' 

2 s w 

m' 
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PT = 1.034 




^2 = -5.745 
F = -12.3 (1-1) 

3^ 



J = -6.37 



«vp 



Z ' 
w 



m' 



The following points are now evident concerning the transfer 
functions defined by equations ( 7 ) and (8): 

1. both gains are always negative since Zw* is a negative 

m' 

number and K ' is always positive, 
s 

2. both are stable transfer functions in that neither has 
a right half plane pole. 

3. (8) is non~ minimum phase since is negative. 

Using these relations, the system to be investigated is re- 
presented by the block diagram of Pig. IV. 




Fig. IV 



(9) 



The open-loop relation for the unity feedback system is; 

(s) 






( s) 



(s) 



K , G T / ^ 

cl cl (s) 



'-^I^c2«l°c2 



TO ^ 2^2 
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where ^c 2 ^c 2 compensation functions to be formu- 

lated, and 



(s) 



KgGg (s) 



(s+W) s 



3 T 2 T T 

s + J_ s + J, 3 + J 
2 1 o 



K, (s'^ + Eg s + P,) 



2 ' “2 ‘ 2' 
(s + N) s 



e(s) = Z* . (a) - Z* (s) 

gi go ^ ' 

Substituting the expressions for (s) and Kg Gg (s) into 

( 9 ), one^Tjfains 



(10) Z’o 



(s) 



^ol ^2°cl ■" V ^2^ 

\ _ - 

s + s + s + J - K, K ^ G ^(s)(s+N)s 
2 1 o 1 c 2 c 2 



Now we choose G _ (s) to be of the form 
c ^ 

2 

G^(s}= s-fas-fli) 

s (a ^ H) 

in which we note that N is a constant and therefore independent of 

•» 

parameter variations. 

With G (s) of this form, (lO) becomes 

C 

(11) ^(s) . Kj (.) 






= * 'V ^l'=c2>= (•^1- >'1^02“)= Wc2'>) 



and the utility of this selection of G o(®) is immediately obvious. By 

c ^ 

suitably choosing K a, and b, the poles of the open-loop relation 

c ^ 

may be conveniently located in the s plane. (As will be indicated 
later, the situation is not quite so straightforward due to saturation 
in 63 .) 



17 



As stated above, the system will be designed with the primary 
parameters assigned their normal or mid-range values. Inserting these 
numbers in (ll) yields 

(12) Z- (s) -233K O (s)(s 2- 5.745s- 11.8) 

go “ 

• e (s) s^ +(5.4 + 2.77K^2^ 2.77K^2^^ ® 

(.382 + 2.77 K^2^^ 

For the present, take equal to unity. Then the zeros 

of the open-loop relations are located at + 7*35 “ 1.6l. Since 

must he positive, we are concerned with variations of the negative 

gain factor -.233 K ^ as far as stability and performance of the closed- 

ol 

loop system are concerned. The locations of the poles of (l2) are 
determined hy experimenting with the root locus plot of this relation 
to achieve the desired resonant frequency and damping ratio of the 
dominant pole-pair of the closed-loop function. The resulting config- 
uration’ is shown in Pig. V. It is to he noted that the most negative 
pole of ( 12 ) is one of the dominant poles. 

With the poles of (12) located at -1, -4, and -9 as shown, the 
denominator of (12) becomes 

s^+ 14s^+ 49s + 36 - 8^+ ( 5.4 + 2.77K-q2)s^ + (3.68 + 2.77 

(.382 + 2o77K^2^) 

Equating coefficients of the various powers of s yields 

'^02 ■ 
a = 5 • 26 

b = 4.14 

The phase-magnitude plot of the open-loop relation incorporating 
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PIG. VII 
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these values is shown in Fig. VI as curve A. Prom this plot it is 

evident that should be a lag network in order to increase the low 

frequency gain of the system. Taking G ^ (s) to be given by 

cl 

G ^ (s) = B + 2 
s + .2 

curve A is modified as shown in curve B of Pig. VI. Applying curve 

B to a Nichol's chart indicates a required gain factor of 14 db for 

an M of 2.3 db. That is 
P 

20 log .233 = 14 

K ' = 21.4 
c 

The frequency response of the closed-loop system is shown in Pig. VII. 

A steady state error of -.4 db, or 4.5%» is indicated in depth rate, 

but this is considered acceptable. 

With K , 0 , (s) and K _ G _ (s) thus determined, the overall 
cl cl c2 c2 ' ’ 

system was simulated on a Reeves Electronis Analog Computer. As will 

be indicated, it was necessary to modify G^^ (s) in order to realize 

a satisfactory step response. The computer ^set-up is described in 

Appendix A, and the plotted responses in § ,0, and e for a step 

s 

input of Z' . are shown in Appendix B, 
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RESULTS 



A. General 

The proposed system was simulated on a Reeves Electronic Analog 
Computer as described in Appendix A. A constraint of plus and minus 
25 degrees or .436 radians was placed on the maximum stern plane angle 
permissible. The input the system was chosen to be a step function. 



Z'i - -5 u.j(r) 



gi 

This choice was made on the basis that in a steady glide w' should be 
desirably zero making 



. 

Z* 

q 



E — 0 



where B is the maximum permissible steady pitch angle. This value of 
0 was taken to be -.5 radian or -28.6°. 30° is usually taken as the 
maximum allowable value of 0 in medium speed submarines. The actual 
peaik value of 0 indicated in the results of Appendix B are only ap- 
proximate due to the assumption that © = sin 0. The error introduced 
by this approximation is Justified in tie sense that using a high 
value of 0 permitted observing the system's general performance in the 
most extreme condition normally encountered. 

B. Evaluation of Results 

With the system as designed in Chapter II set up on the computer, 
and the primary parameters assigned their normal values, it was observed 
that the overshoot in 0 and e were of the order of thirty percent. This 
was due to excessive lag introduced by and saturation in ^ Ac- 

cordingly, the locations of the pole-zero pair of the lag network were 
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shifted to -0.1 and -1.0 respectively by changing the setting of po- 
tentiometer number I 5 . This corresponded to changing the value of 
k in the lag function representation given in Appendix A from 1 to .5* 
The effect of this adjustment on the phase-magnitude plot of the open- 
loop function and on the closed-loop response are shown in Curve C of 
Pigs. VI and VII. After this adjustment, the system performed ’satis- 
factorily. Plotted results are shown in Appendix B. 

The following observations are to be noted concerning these 
re suits: 

1. the system specifications are met for values of Z^, of 

m' 

-I .5 and -2.0 for all values of U greater than 2.0 and 
all values of considered. 

2. the system does not function satisfactorily for Z ’ =-1.0. 

w 

m* 

At values of U of loO and 2.0, the system either responds very 
slowly and with large steady-state error, or shows very little res- 
ponse at all, especially for low values of K*. This is due directly 
to the small lift force generated by the stern planes at slow speed 
and cannot be taken as a failure of the plane angle control system. 

At these speeds, S takes on a maximum value and retains this value 
as would be expected. For this reason values of U less than 3.0 are 
not considered in the results. 
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IV 



DISCUSSION OP RESULTS 

It was noted that the system did not function satisfactorily 

for 2^/®' equal to -1.0. The reason for this becomes apparent 

when the effects of this parameter and of saturation is <5^ on the 

locations of the open-loop poles of the system are considered. 

Limiting the maximum value of 5 is equivalent to limiting the ef- 

s 

fectiveness of the compensation gains and is recalled 

that the open-loop poles of the system were moved further out into 

the left half of the s p].ane by a suitable choice of effect 

of saturation in 5 is, therefore, to constrain these poles to lie to 

s 

the right of their desired locations. The effect is aggravated by 

the tendency of as well as to approach zero as U increases when 

Z '/m* equals -1.0. For values of Z '/m' between zero and -1.0, J- 

approaches a negative number. In the extreme case, two of the poles 

fall between the origin and the zero of the lag network, and these 

become the dominant pole pair. The net result*’is a lowering of the 

closed-loop natural frequency and damping for a given gain. Two 

possible remedies present themselves. First, the pre-saturation gain 

may be decreased to keep the system in the linear region of 5^ 

longer, with a resultant decrease in steady-state accuracy. Pig. XXVI 

of Appendix B shows this result for 9»25 instead of 21.4. A 

second approach lies in increasing the post-saturation gain of the 

hull via K* . This is equivalent to supplying more control surface 
s 

area. Pigs. VIII through XIII of Appendix B indicate this tendency. 
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A more detailed discussion of the effects of saturation indicated 
would require a non-linear analysis of the system and will not he 
attempted here. 

The question immediately arises as to why these effects were 
not predominant at the higher values of investigated. The 
necessity of shifting the zero of the lag network to the right men- 
tioned above indicates that these effects were present when Z '/m’ 
equalled -1.5> ^ut a suitable remedy was quite simply achieved. The 
effect is not so pronounced simply because and assume larger 
values and the hull gain is larger for these values of Z^'/m'. 

This failure of the system at low values of Z ' is not considered 

w 

a serious limitation on the basic approach presented. For normal hull 

/ 15 

forms, the magnitude of Z^,/m' is rarely less than 1.2. The re- 

sults simply confirm the known fact that hulls with little directional 
stability must be provided with larger control surfaces. 

The remainder of the plotted results of Appendix B for values 
of Z^'/m' of -1.5 -2.0 are considered to justify the method pre- 

sented. The behavior of pitch angle, a critical quantity, in the dive 
entry maneuver simulated is considered to be quite satisfactory. The 
error in steady state depth rate is maintained less than eight percent 
in all instances, and the overshoot is negligible. It will be re- 
called that the calculated steady state error of the system as com- 
pensated was four and one-half percent. Prom Pig. XVII, the steady 
state error of the system as designed is seen to be four percent. Since 
U appears in the denominator of the expressions for and J^, decreas- 

a/ 

ing U will move the poles of the open-loop relations to the left and 
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thus increase the steady state error, an effect clearly shown in the 
plotted results. The converse is true when U is increased. It is 
seen from the plots that this effect is made quite small by the 
compensation technique employed. 

It will be recalled that the requirement for a steady-state 
plane angle to counteract the metacentric righting moment was neg- 
lected in the system formulation. Prom the results, it is seen that 
this requirement is almost completely satisfied by the system due to 
the effect of metacentric righting moment on depth rate. For low 

r>J 

values of U, an appreciable steady state plane angle exists but 
there is still a small pitch rate. 

It is noted that the effect on steady state error due to changes 
in hull gain through variations in K* and Z smaller than the 

s w ' 

effect of metacentric moment noted above. However, several runs were 
made with adjusted to correct for this change in hull gain. The 
runs were made only for the case U = the results shown in 

Pigs. XXVII and XXVIII indicate that a slight improvement in steady 
state error results. 

In the basic system presented, it was assumed that an infinite 

plane rate,? , was available. This is not the case in practice. The 
s 

compensation method presented, however, can easily be extended to 

handle the case when plane actuator dynamics are included in the 

function. For example, if these dynamics are representable by a 

transfer function of the form ^ , addition of the term (p s+l) 

p s+l 

in the numerator of both G , and G „ would cancel this function and 

cl c2 

require minor changes in the constants associated with 
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^c2^c2* -A^gain, complete cancellation would be achieved only in the 
linear region of stern plane operation. 

It is becoming common practice to evaluate performance of a 
proposed submarine control system in terms of the comparison between 
the trajectories achieved for step inputs to the system and experi- 
mentally determined minimum time or limit trajectories. Due to time 
and computer capacity limitations, such a direct comparison was not 
carried out in this investigation. It is maintained, however, that 
on the basis of the indicated behavior of stern plane angle the 
plotted trajectories of depth rate error and pitch angle for the 
cases of interest closely approach those observed in the limit dive 
entry. In the results of Appendix B for the cases in which Z^*/m* 

assumes the values of -1.5 and -2.0, § shows at most one complete 

# s 

reversal from its initial maximum value "before it quickly settles 
out to its steady state value. This behavior is necessarily charac- 
teristic of a minimum time trajectory. More will be said on this 
point in the following chapter. 
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V 



A PROPOSED PULL-OUT TIME PREDICTOR 

The proposed system functions satisfactorily in carrying the ship 

through a dive or rise entry phase and maintaining the required dive 

or rise, hut it cannot inform a human operator when to pull out to 

achieve a desired depth. As was indicated, a computer or programmer 

could be substituted for the human operator, but a simpler solution 

would be desirable. In the course of the investigation reported 

herein, a simple analog method of computing pull-out time presented 

itself. The basic idea stems from an optimum controller concept 

proposed by Lieutenant M. C. Davis USW in his doctoral thesis on multi- 

17 

dimensional random processes. This technique is considered worthy 
of note and further investigation as an auxiliary component of the 
proposed system. 

If the optimum trajectory in depth is specified to be one executed 
in minimum time with no overshoot, then the optimum controller applies 
maximum control effort throughout the entry and pull-out phases. 
Consider a ship controlled by the proposed system to be in some steady 
dive. In pull-out, the operator could simply return his control col- 
umn to the zero position, commanding zero depth rate, or he could com- 
mand varying amounts of negative depth rate to speed up the pull-out. 
But we have seen that the system presented in this paper approaches 
the performance of the optimum controller itself when the input is the 
desired step change in depth rate. Therefore, commanding a negative 
depth rate in pull-out cannot give a faster response than simply 
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returning the control column to zero. The problem, then, is deter- 
mining when the operator should zero his control column to level out 
at a desired depth. 

Th problem and an indication of a possible solution are shown in 
the following figures. In (a), the ship is shown in some steddy dive. 




ib) 



as indicated by the solid line, which will carry it through the de- 
sired depth At successive instants, t^, the present state of 

the system cas determined by speed, depth, depth rate, pitch, and 
pitch rate are introduced as initial conditions into a fast time 
analog of the system, shown in (b). The input to the analog is a 
step function of equal magnitude but opposite sign to the existing 
commanded depth rate, a step function returning commanded depth 
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rate to zero. The computing periods, in real time, are short relative 
to the response time of the hull. The output of the analog is inte- 
grated to yield, in fast time, the trajectory in depth which would 
result if the control column were zeroed at the associated sampling 
instant t^. This projected depth response is subtracted from the 
desired depth, represented by a constant quantity, and the difference 
is supplied to a zero-crossing detector. At the first zero-crossing 
of this difference, the indicator is actuated. Now it is noted from 
the dashed representations of the projected trajectories in (a) that 
at the first zero in this difference both the depth error and its 
first derivative are zero. That this will be the case is indicated 
by the plots of depth rate error response in Appendix B. This simp- 
lifies the operator's task since if he zeroes his control column at 
the first indication of a zero crossing he will come out on depth 
with depth rate passing through zero. Then he need only correct for 
the error in depth due to the small overshoot of the depth rate con- 
trol system. Note that pullout time is based on the projected depth 
at peak overshoot in depth rate so that in a dive this error will 
always be on the safe side. 
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VI 



CONCLUSIONS AND RECOMMENDATIONS 

On the basis of the results shown in Appendix B, it is concluded 
that the system presented is entirely feasible and worthy of further 
investigation and implementation. It is noted that this system pos- 
sesses the following very desirable characteristics from the operator's 
viewpoint: 

1. the system adjusts itself to any steady dive, or rise 
commanded by the operator in a minimum time. 

2. response becomes faster as speed increases, but does 
not become oscillatory or unstable. 

3. pitch angle takes care of itself within the system, 
leaving the operator concerned only with the primary 
quantity, depth rate. 

4. the control system minimizes the effects of speed and 
hull geometry variations by increasing the natural 
frequency of the control loop by a factor of at least 
five. 

It is considered that this system, in conjunction with a pull-out 
time indicator discussed in the proceeding chapiter, will provide near 
optimiim depth changing and depth-keeping capability. Due to its 
broadband nature, one system of this type will provide adequate con- 
trol characteristics in a variety of hulls, obviating the need to 
tailor-design a control system for each new ship. 
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APPEITDIX A 



COMPUTER STUDY 

A. Computer Simulation of the Proposed System 

The system as shown in Fig. IV was simulated on the analog 
computer in the following manner. To avoid loosing sight of the 
basic problem, the force and moment equations of the hydroframe as 
given in equations (3) and ( 4 ) were used to represent the hull. To 
facilitate the computer solution, these equations were rewritten in 
the following form 



(3a) 100 



e ■ f 



100 - 10 . 1 1 z. 

dr 



+ .1(100^') 



( 4 a) 100_dfo_= .123 |;^qq^, _ 1000 - 



' 1 f.5 (100 dO) + K' (10 ) + 

I I s s 

' I dr 

']] 



where 



r« 



dr 



10'^ 






u 



10 .11 w 

1 m' 



.25 (1000 1 ^) + 



+ .5 (lOOw') + .5 K' (10)(10 ) 

s s 



']} 



d 

d 2 r 



10 



dr 



100 



dt" 



dr 



A point concerning the analog treatment of the lag network is 
worthy of note. Writing the defining relations of this component in 
the following form. 



.?Bl 

e ( s) s + .2K 



and carrying through the indicated operations and change of time base 
yields 

10 Ifli = 1C (10 1^) . K (2K . - .2 ) 

dr 
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is unavailable, and noting that the 



Integrating once since 4^ 

dr 

initial conditions on the integration are zero yields 

*31 “ ■'ol (1° + kJ" e - .2 SJ dr 

or - ‘Ol^cl 

Now the coefficient of the integral will appear in the computer 
diagram simply as a potentiometer for values of K between 0 and 10. 

K determines the location of the pole-zero pair of the lag network 
along the negative real axis of the s plane relative to the selected 
position at -2 and -.2, the ratio of their distances from the origin 
remaining constant. This facilitates a simple method of experimen- 
tally varying the phase shift introduced by the lag network without 
affecting the steady state error of the system. 

The feedback compensation is represented by the following; 



10 ^ 

dr 



.310 (1000 ) + 

2 



dr‘ 



.163 (1000 ) + .128 (100 e) - 

dr 



-.104 (10 6^2 ) 

A limiter is employed to restrict the maximum value of S g 

to + 25 ° or .436 radian. 

Z' is computed from the relation 
go 

o 

100 z' = 100 w’ - 100 e 
go 

The complete computer diagram is shown on Plate I. 



31 




RERC CIRCUIT 



APPENDIX B 



ANALOG COMPUTER RESULTS 



Pigs. VIII through XXV show the results of the analog computer 

study of the proposed system as plots of tf , &, and e vs the non- 

s 

dimensional time parameter for various values of the three primary 



parameters K^, 


, and 


U. The values of these parameters con- 


sidered are as follows: 




-V 

U 


3.0, 


5.0, 7.0, 9.0 


K * 
s 


C\J 

• 

o 


0.3, 0.4 


Z • 






w 


-1.0, 


-1.5, -2.0 


m’ 






Curves for each value of 


U are shown on one plot for each pair of 


the parameters 


and 


Z ' considered, 

w 



m' 



Pig. XXVI shows the performance of the system in pitch angle 

and depth rate for Z ’/®* “ -1«0, U = 5>0, K ' = 0.3, and K , =9.25 

w' ’-''s 'cl 

Pigs. XXVII and XXVIII show this data for Z * of -1.5 and -2.0 



w 



m' 



and all values of K^' considered with adjusted to maintain the 



product constant. 
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